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Abstract 



We adopt a personal approach here reviewing several calculations over the 
CN ■ 

^ ■ years in which we have experienced confrontations between cluster models and 

the shell model. In previous cluster conferences we have noted that cluster 
models go hand in hand with Skyrme Hartee-Fock calculations in describing 



■ states which cannot easily, if at all, be handled by the shell model. These 

I are the highly deformed (many particle - many hole) intruder states, linear 

chain states e.t.c. In the present work we will consider several topics; the 
Q ■ quadrupole moment of ^Li, the non-existence of low lying intruders in ^Be, and 

pH I then jumping to the /7/2 shell, we discuss the two-faceted nature of the nuclei 

- sometimes displaying shell model properties, other times cluster properties. 

X 

^ ■ 

^ I. THE QUADRUPOLE MOMENT OF THE J = 1+ STATE IN ^LI 

Whereas the quadrupole moment of the deuteron is positive(Q = +2.74mb), that of the 
J=l+ state of ^Li is negative, Q=-0.818(17)mb. The magnetic moment of the deuteron is 
/i = 0.85741 nm while that of "^Li is 0.822 nm. 

There appears to be a big discrepancy between cluster model calculations and the shell 
model calculations. In nearly all cluster model calculations Q comes out positive. However 
in many shell model calculations Q comes out negative, sometimes too negative. This is 
an important problem that deserves further attention. See for example arguments in the 
literature between the cluster group [1] and the shell model group [2]. See also the recent 
compendium of A=6 by D.R. Tilley et. al.[3]. 
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For example in a modern shell model approach by Forest et. al. [4] gets about -8 mb 
for Q, a factor of 10 too large but of the correct sign. On the other hand in a dynamical 
microscopic three cluster description of ^Li where the clusters are a, n, and p the result is 
Q = 2.56 mb.[l] 

In shell model calculations that we performed [5] we started with 2 valence particles in 
the Op shell (Ofiw). Then we allowed up to 2 ?iw and then up to 4 hu; excitations. In the 
hcu space if you do not have a tensor interaction Q comes out positive. With a 'realistic' 
tensor interaction Q comes out negative but too negative Q=-3.5mb. However with a former 
student Zheng, who at Arizona also developed the no core approximation with Barrett et. 
al.[6], we showed that when higher shell admixtures were admitted Q became smaller in 
magnitude and closer to experiment as shown in the following table. [5] 

The results are shown in the following table 



SPACE 


Q(mb) //(mm) 


hu 


-3.60 0.866 


2 hw 


-2.51 0.848 


4 ?iu; 


-0.085 0.846 


Experiment 


-0.82 0.822 



Note that the shell model calculations cannot get the magnetic moment low enough. 
With up to 4 huj admixtures wc actually overshoot and get a quadrupole moment that 
is too small but still negative. Some cluster models appear to explain the low magnetic 
moment. 

An excellent discussion of many shell model calcultion of Q and n has been given by 
Karataghdis et. al [7]. The value of Q that they obtain with what they call the "Zheng" 
interaction [8] in the up to 0,2,4 and 6 tiu; spaces are -2.64,-2.08,-0.12, and 0.17 mb respec- 
tively. Thus they get Q to become positive at the 6 fiu; level. But then they quote Zheng et. 
al [8] as getting a value of -0.67mb in the same 6 fiw space. It is not clear why the two cal- 
culations give different answers. The changes in // in ref [7] are more moderate 0.869,0.848, 
0.845 and 0.840 nm in the up to 0,2,4 and 6 fiu; spaces. 

Looking at all the calculations by all groups (including our own) , the situation is certainly 
confusing, and the problem deserves further attention. This is certainly a basic problem, the 
deuteron embedded in the nuclear medium. This problem has wider implications whether 
or not there is T=0 pairing can depend on how higher order configurations affect the tensor 
interaction in the valence space. 
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II. ABSENCE OF LOW LYING INTRUDERS IN «BE AND THE a PARTICLE 

MODEL 

The 0"*" bandhead for low lying intruders in ^^0, ^^C, and ^°Be are at 6.05 MeV, 7.65 
MeV and 6.18 MeV respectively. In ^^0 and ^^C, these are predominantly 4 particle 4 hole 
excitations. In ^^C, we identify the intrinsic state as a linear chain. In the 7th edition of the 
'Tables of Isotopes' possible intruders in ^ Be were indicated, a J^O"^ state at 6 MeV and 
J=2+ state at 9 MeV. 

In shell model calculations allowing 2 particle 2 hole excitations we were able with a 
quadrupole-quadrupole force to get a J=0+ state at 9.7 MeV in ^°Be, too high but in the 
right ballpark. [7,8] But we could not get low lying intruders in ^Be below 30 Mev.[7,8] We 
used a deformed oscillator model to show why one gets intruders in ^^C and ^°Be but not ^ 
Be. 

But perhaps the simplest explanation as suggested to us by E. Vogt is given by the a 
particle model. In ^^C we can rearrange the a particles from a triangle to a linear chain. In 
^Be we have only 2 a particles. One can get a rotational band by having the 2 a's rotate 
around each other but that is all. 

The mere existence of these intruder states is of astrophysical importance. In the beta 
decay ^B— >^Be + e"*" + z/ one goes from a J=l+ T=l to J=2"'" states. This is the famous 
'Ray Davis' neutrino. If there were a 2+ state at 9 MeV then there would be more high 
energy a's than there would be if the decay were to the 2i state at 3.04 MeV. [9] The alpha 
spectrum from the decay of ^Be seems to show more high energy alphas, but we would say 
that they are not due to low lying intruders. 

III. CLUSTERING AND SHELL MODEL IN THE Fy2 REGION 

In a previous cluster conference in Santorini (1993) a spectrum of ^^Ti was shown in 
an a cluster model. [10] The spectrum looked reasonable except that there was a wide gap 
between the 10"*" and 12+ states. However, these states arc sufficiently close together that 
the 12+ state is isomeric. In a single j shell basis (j=f7/2) ^^Fe is the 4 hole system and 
it should have an identical spectrum to that of ^^Ti provided the same interaction is used. 
However in ^^Fe the 12+ hes below the 10+ state. It is extremely isomeric and has a hfetime 
of 12 minutes. 

We have studied this and other topics by calculating the spectrum of ^Ti (^^Fe) with a 
variety of interactions designated as Model X. (See Tables I, II) 
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Model I: Use the spectrum of ^^Sc as input (particle-particle) Identify < {j'^)'^V{j'^y >— 
E{J) experimental. For isospin T=0 J can be 1,3,5 and 7 while for T=l J is even 0,2,4, and 
6. 

Model II: Use the spectrum of ^^Co as input (hole - hole). If there were no configuration 
mixing these two spectra would be identical. However, there are some differences eg the 1^ 
state is much lower in ^^Co than in ^^Sc. 

Model III: Now we play games. We want to find out how important are the T=0 matrix 
elements for the structure of the nuclei, (e.g. Is T=0 pairing important?) Noticing that in 
^^Sc the J=2,3 and 5 states are nearly degenerate in this model we set all the T=0 matrix 
elements to be the same and to all equal -E(2+) = 1.5863 MeV. 

In model III we then have V'^=^ = V{^'^Scf=^ J =0,2,4,6 and y^=° = constant^ E{2+) 
J=l,3,5,7. We can then write \^'^^°=constant (1/4 — ti ■ ^2)- We can then write V'^^^ = 
c(l/4 — ti ■t2) where c is a constant. Hence X^j^^- y^J"'^ = c/8(n(n — 1) + 6) — c/2T(T + l). This 
means that the spectrum of states of a given isospin e.g. T=0 in "^^Ti(^^Fe) is independent 
of what the constant is, it might as well be zero . Of course the relative splitting of T=l and 
T=0 states will be affected. Model HI will be the standard from which we derive Model IV. 

Model IV: Relative to the degenerate case above, we now move the J = 1^ state down 
in energy to 0.5863 McV. Our motivation is based on numerous discussions about the im- 
portance of T=0 S=l "pairing" in nuclei. We hope to simulate the T=0 pairing by this 
lowering. 

Model V: Relative to Model HI we bring the J=l"'" and J^?"*" states down to an energy 
of 0.5863 MeV but keep the J=3+ and 5+ at E=E(2+)=1.5863 MeV. This spectrum is very 
close to that of ^^Sc. 

IV. DISCUSSION OF RESULTS 

Let us first compare Model HI (all T=0 matrix elements are degenerate) with Model 
I (spectra of ^^Sc). As already mentioned, making T=0 matrix elements degenerate is 
equivalent to making them zero as far as T=0 states are concerned. 

The main difference is that the states with J=6,4,7, and 8 come down in energy as does 
J=9+. Also the 12-10 gap is a bit greater than for the ^^Sc spectra case, reminiscent of 
the a particle model. The J^Q"*" state is below the lO"*" and 12+ in the degenerate case. 
Clearly it is the high energy side of the spectrum which is most sensitive to the change from 
experimental spectrum to the "T=0 degenerate" case. 
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Despite the changes, we can say that the T=l two body matrix elements give the domi- 
nant structure of the spectrum whilst the T=0 matrix elements provide the fine tuning. 

We now compare Model IV with Model III. The only difference is that we break the T=0 
degeneracy by lowering the J^l"*" state from 1.5863 MeV to 0.5863 Mev. We hope that this 
simulates to some extent T=0 S=l neutron-proton pairing. The change from degenerate 
case is not that large. There is a tendency to go towards the spectrum of Sc. The J=3,5,7, 
and 8 states are raised somewhat in energy. However it is hard to find a clear signature of 
this S=l pairing. 

Not shown is Model V where we bring down both the J=l+ and 7+ states to 0.5863 MeV 
keeping J=3 and 5 at E(2) = 1.5863 MeV. This input spectrum is close to that of Sc so 
that it is not surprising that the ^^Ti spectrum is likewise close. 

We lastly consider the results using the spectrum of ■'''^Co. This was done some time ago 
by Geesaman [11,12] Note that there are significant changes, all at the high energy high 
angular momentum part of the spectrum. Relative to the "^^Sc case the 10+ and 12"*" states 
are down in energy with the 12+ below the 10+ thus leading to a long lifetime for the 12+ 
state. Note that the 9+ state is now at a much high energy than the 10+ or 12+. Recently 
the 10+ state in ^^Fe, which lies above the 12+ has been found in ^^Fe by Ur et. al. [13]. It 
would also be of interest to find the 9+ state. 

V. MANY PARTICLE, MANY HOLE STATES IN ^°CA 

This is a topic we discussed in previous cluster meetings so we will be brief. [14] We 
just want to remind the reader that there are all sorts of many particle-many hole highly 
deformed states in ^°Ca. One cannot properly describe ^°Ca in a cluster model consisting 
of ^^Ar plus an alpha particle. At the very least one has to start with ^^S plus two alpha 
particles. 

In a Skyrme Hartree-Fock calculations (SK III) we obtain a near degeneracy of the 4p- 
4h and 8p - 8h intrinsic state. The respective energies are 12.1 MeV and 11.4 MeV. The 
8p-8h intrinsic state energy is lower than the 4p-4h. By the time projection and pairing are 
included, the 4p-4h comes lower than the 8p-8h (6.85 MeV vs 8.02 MeV) in agreement with 
the order the J = excitation energies of 3.0 and 5.1 Mev. Pairing will lower the states 
even more. We actually found many more deformed states of the form np-nh n=2,3,4,5,6,7 
and 8. The intrinsic states are nearly degenerate in energy - we called this a deformation 
condensate. We also found for these states that the deformation parameter was approxi- 
mately proportional to n i.e. the value of (3 for 8p-8h is approximately twice the value of (5 
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for 4p-4h. 

In ^°Zr one of the np-nh states becomes the ground state. This is the 12p-12h state 
which has a calculated value oi P — 0.4 A more superdeformed 16p-16h state with (3 — 0.6 
is calculated to be at an excitation energy of about 8 MeV. 

VI. TWO DIFFERENT VIEWS OF THE Fy2 REGION 

In March 2000 issue of the Physical Review C 61 there are two papers side by side. One 
is by our group [15] and one by H. Hasegawa and K. Koneko [16]. We both do calculations 
in the f7/2 shell. We emphasize shell model behavior whilst the other authors the a cluster 
behaviors, even though their model space is limited to f7/2. 

The other authors point out that we can get an excellent approximation to the ground 
states of Up = n„=2m nuclei (n^ is the number of protons e.t.c). 

|(/7/2)'"^/ = T = >= -^{alr\Ao > (1) 
where (aj) creates a 4 nucleon cluster. 

al = 5:(Jt, Jr : / = T = 0)(4,4J,=ot=o (2) 

J,T 

For '^^Cr this approximation give -32.04 MeV for the ground state energy where as the exact 
value is 32.70 MeV. 

We on the other hand have emphasized the shell model aspects. [15] In the previously 
mentioned paper, we find an approximation for the excitation energies of single and double 
analog states in the f7/2 region and in a an earlier paper "Fermionic Symmetries: Extension 
of the two to one relationship between spectra of even even and neighboring odd mass 
nuclei" [17] we noted two things. 

A. There is often a two to one relation between spectra of even-even and even odd nuclei, 
and in some cases the single j shell model predicts this. 

B. Excitation energies of analog state are approximately the same if the neutron excess (or 
equivalently the ground state isospin) is the same. 

The above results can be parametrized by the following formulae 
SINGLE ANALOG EXCITATION (SA) 

E{SA) = b{T + X) (3) 

DOUBLE ANALOG EXCITATION (DA) 
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E{DA) ^2b{T + X + l/2) 



(4) 



This formula will give a two to one ratio for E(DA)/E(SA) for (^^Ti,^^^^)^ {^^Cy,^^Cy), 

The experimental SA and DA are shown in the following table. 

T=0 4^Ti (9.340), ^^Cr (8.75), ^^Fc (8.559) 

T=l/2 ^^Sc (4.274)'', ^^Ti (4.338)", ^^Ti (4.176), 49Cr(4.49), 

^^Mn (4.451), ^^Co (4.390), ^^Fe (4.250) 
T=l ^6rp- (14.153)^ 50^^ (13.222) 

T=3/2 ^^Sc (6.752)^ ^^Ti (7.187), ^^Cr (6.611) 

T=2 4^Ti (17.379) 

T=5/2 4^Sc (8.487)", ^^Ti (8.724) 

a - obtained from binding energy data. 

In Tabic II we compare the theoretical single j shell calculations with the linear formula. 
We take b = 2.32 MeV X=1.30. Note that in the SU(4) limit X=2.5. The fact that 
SU(3) gives the hnear formula is not sufficient for it to be the correct theory. For a simple 
monopole-monopole interaction a+bt(l)t(2) X=l. 

Some of the two to one ratio's hold rigoursly in the single j shell model. This holds for 3 
particle and 4 particle systems or 3 holes and 4 holes, eg (43Ti,^Ti), (^3Sc,^^Ti),(53Fe,52pg) 
Here not only single or double analog but all the J=j states in the odd spectrum are at half 
the energy of the J=0+ states in the even system. 

Some of the relations hold approximately in the single j shell model eg for (^^Sc,^^Ti) and 
for the cross conjugate pair (^^Cr,^°Cr) we would get a 2 to 1 ratio if seniority four states 
could be neglected and one only had v=0 and v=2. 

Miraculously the 2 to 1 ratio holds remarkably well experimentally for (^^Cr,^°Cr) - the 
values are 6.511 and 13.022 MeV respectively, this despite the fact that in the single j shell it 
should only hold approximately. Ironically the simplest system for which 2 to 1 should hold 
exactly does not work so well. That is to say for (^^Ti,^'^Ti) the values arc 4.338 and 9.340 
MeV. When configuration mixing is included, agreement with the deviation is explained. 
This might be an example of a 4 particle clustering. For the hole system (^^Fe, ^^Fe) on the 
other hand 2 to 1 works much better. 

The fact that there is in general a close relation between even-even and even-odd puts 
to question whether in those many cases there is any a particle clustering. 

The single j shell calculation does not predict exact equality for the S.A. excitation 
energies in "^^Sc and "^^Ti. The relative values are very close however 4.142 and 4.112 MeV 
respectively. This is fascinating. We take ^^Sc, jam a deuteron into it to form ^^Ti and it 
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seems hardly to make any difference for S.A. excitations. [15] 

VII. TWO VIEWS OF CROSS CONJUGATE RELATIONS 

In the single j shell model, the spectra of cross conjugate nuclei should be identical. [18] 
(for j" states.). A cross conjugate nucleus is one obtained by changing protons into neutron 
holes and neutrons into proton holes. The cross conjugate of ^^Ti is ^°Cr. Let us compare 
the spectra 



0.000 0.000 

2 0.889 0.787 0.8853 

4 2.010 1.884 0.9373 

6 3.297 3.164 0.9597 

8 4.896 4.740 0.9681 



The fits are very good. The ^'^Cr excitations are slightly smaller - it could be a universal 
A dependence. Where does the remarkable agreement leave room for a clustering? 

However we can look for other things besides the spectra. In a recent experiment theory 
collaboration where the leading experimentahsts were N. KoUer and H.A. Speidel [19] good 
agreement was obtained for g(2+) in ^°Cr but bad agreement for ^^Ti. The shell model 
predicts a high value of g(2+) and g(4''") for ^°Cr but low values for ^^Ti 0.25 mm. The 
high values are confirmed for ^°Cr but for ^^Ti the measured g factors are closer to 0.5 which 
suggest the rotational value gR—Z/ A. These results suggest that there must be considerable 
clustering in ^^Ti that is not present in ^°Cr. In general the shell model appears to work 
better in the upper half of the "f7/2 shell" than in the lower half. There appears to be much 
more going on in the lower half and probably this is due to intruder /cluster mixing in with 
the basic shell model states. 

VIII. CLOSING REMARKS 

We have provided several examples where Cluster Models and the Shell Model confront 
each other usually to the mutual benefit of both models even though in the short term 
there might be some arguments. The two models give the opposite sign for the quadrupole 
moment of ^Li, and this has to be resolved. The cluster model provides insight into some 
results of detailed shell model calculations e.g. why there are no low lying intruders in 
^Be. The low lying intruder states e.g. 4p-4h and 8p-8h in "^^Ca are essentially impossible 
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to calculate in the shell model. However here cluster models and Skyrme-Hartree Fock go 
together in describing such states. In the £7/2 region we raise the question (without fully 
answering it) of how to distinguish symmetry energy from clustering energy. Finally we 
point out the issue of hidden clustering. 

This work was supported in part by the U.S. Department of Energy DE-FG02-95ER- 
40940 
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TABLES 



TABLE 1. Two particle Matrix Elements < if yV{fy > 



Case 


.1=0 T=l 


J=2 T=l 


J= l T=l 


J=6 T=l 


.1=1 T=0 


.1=3 T=0 


.1=5 T=0 


,1=7 T=0 


Model P) 


0.000 


1.5863 


2.8153 


3.2420 


0.6110 


1.4904 


1.5101 


0.6163 


Model 11^) 


0.000 


1.4465 


2.6450 


2.9000 


0.9372 


1.8224 


2.1490 


0.1990 


Model IIP) 


0.000 


1.5863 


2.8153 


3.2420 


1.5863 


1.5863 


1.5863 


1.5863 


Model IV'') 


0.000 


1.5863 


2.8153 


3.2420 


0.5863 


1.5863 


1.5863 


1.5863 



TABLE IL ^^Ti {^'^Fe) spectra for Model I,II,III and IV 



Model I (J'^ E)") 


Model II (J'^ E)^) 


Model III (J'^ E)=) 


Model IV (J'^ E)'^) 


0+ 


0.000 


0+ 


0.000 


0+ 


0.000 


0+ 


0.000 


2+ 


1.163 


2+ 


1.015 


2+ 


1.303 


2+ 


1.253 


4+ 


2.790 


4+ 


2.628 


4+ 


2.741 


4+ 


2.800 


6+ 


4.062 


6+ 


4.079 


6+ 


3.500 


6+ 


3.738 


3+ 


5.786 


8+ 


5.772 


3+ 


4.716 


3+ 


5.031 


5+ 


5.871 


7+ 


6.018 


5+ 


4.998 


5+ 


5.082 


7+ 


6.043 


12+ 


6.514 


7+ 


5.356 


7+ 


5.687 


8+ 


6.084 


3+ 


6.540 


8+ 


5.656 


8+ 


6.045 


10+ 


7.384 


5+ 


6.602 


9+ 


7.200 


9+ 


7.731 


12+ 


7.702 


10+ 


6.722 


10+ 


7.200 


10+ 


7.731 


9+ 


7.984 


9+ 


8.048 


12+ 


7.840 


12+ 


8.371 



a) Input is spectrum of ^^Sc(particle-particle) 

b) Input is spectrum of ^^Co(hole-hole) 

c) The T=l matrix elements are from the spectrum of ^^Sc. The T=0 matrix elements are 
degenerate at 1.5863 MeV. 

d) Same as model 3 except that the J=1+T=0 energy is lowered to 0.5863 MeV. 
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